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Abstract 


A lead scintillator sandwich sampling calorimeter has been installed 
in the HERA tunnel 105.6 m from the central ZEUS detector in the 
proton beam direction. It is designed to measure the energy and 
scattering angle of neutrons produced in charge exchange ep colli¬ 
sions. Before installation the calorimeter was tested and calibrated 
in the H6 beam at CERN where 120 GeV electrons, muons, pions 
and protons were made incident on the calorimeter. In addition, the 
spectrum of fast neutrons from charge exchange proton-lucite colli¬ 
sions was measured. The design and construction of the calorimeter 
is described, and the results of the CERN test reported. Special at¬ 
tention is paid to the measurement of shower position, shower width, 
and the separation of electromagnetic showers from hadronic showers. 
The overall energy scale as determined from the energy spectrum of 
charge exchange neutrons is compared to that obtained from direct 
beam hadrons. 


1 Introduction 

The ZEUS[2, ^ collaboration has installed a Forward Nentron Calorimeter 
(FNC) in the HERA tnnnel 105.6 m from the central detector in the down¬ 
stream proton direction. The calorimeter, which views the interaction point 
from zero degrees, is designed to stndy nncleon charge exchange reactions 
in ep collisions by detecting energetic nentrons prodnced at small scattering 
angles. The feasibility of the idea was tested snccessfnlly with two small test 
calorimeters in 1993]^ and 1994[Q. High energy nentrons were detected in 
both photoprodnction and deep inelastic scattering]^ H, |[. 

Here we report on the design and constrnction of the FNC. We present 
resnlts from the test and calibration stndies at CERN, with special attention 
paid to the measnrement of shower position, shower width, and the separation 
of electromagnetic showers from hadronic showers. At HERA the overall 
energy scale is determined from the energy spectrum of neutrons produced 
by charge exchange of beam protons with residual gas in the beam pipe. A 
test of this procedure was made at CERN by comparing the neutron energy 
spectrum from the charge exchange of beam protons in a Incite target with 
that obtained from the direct beam. 
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2 Design and Construction of the FNC 


2.1 Introduction 

The FNC detects neutrons produced at a scattering angle less than a millira- 
dian and with energies as high as the HERA proton beam energy of 820 GeV. 
In physics studies of nucleon charge exchange, the energy transferred from 
the incoming proton [AE = Ej, — E, the difference between the beam en¬ 
ergy and the scattered neutron’s energy) is an important quantity. Excellent 
resolution on E is needed to obtain good resolution on AE when E is near 
the kinematical end point Eh- The best energy resolution for high energy 
hadrons is obtained by the use of a ‘compensating’ calorimeter whose re¬ 
sponse to energy deposited by electromagnetic showers equals its response 
to energy deposited by hadronic showers {e/h = 1). The other physics re¬ 
quirement is a measurement of the transverse momentum, of the neutron. 
This can be achieved by transverse segmentation of the calorimeter and its 
readout. 

To minimize costs, as well as development and construction time, we 
chose to design the FNC as a lead-scintillator sandwich calorimeter read out 
by wavelength shifting light guides and photomultiplier tubes. A calorime¬ 
ter of similar design built by the ZEUS calorimeter group was shown to 
compensate 0. That calorimeter had a unit cell consisting of 1.0 cm of lead 
and 0.26 cm of scintillator. In hadron beams up to 100 GeV at CERN it 
proved to have an energy resolution of 6E/E = 0.45/a/E with a very small 
constant term. After correcting for transverse leakage of hadronic energy it 
had an electron to hadron ratio e/h close to unity. 

The containment of 95% of the energy of an 800 GeV hadronic shower 
requires 9-10 interaction lengths of absorber^, Q. To ensure this depth in the 
FNC we increased the ratio of lead to scintillator thickness from 4:1, which 
gives a compensating calorimeter, to 5:1, in order to meet the space restric¬ 
tions in the HERA tunnel (see section |2.2|) . As a result, we expect the FNC 
to overcompensate {e/h < 1), at least for energies below 100 GeV; however, 
Monte Garlo studies[|T3 with FLUKA indicate that the e/h ratio depends 
only weakly on the sampling fraction. Furthermore, FLUKA predicts a weak 
dependence of e/h on the neutron energy. Hence the FNG should be close 
to compensating throughout the energy range up to 820 GeV. 

Since the sampling term in the energy resolution function increases as the 
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square root of the sampling fraction we expect a degradation of the energy 
resolution by 12%. More serious is the degradation due to transverse energy 
leakage. The Monte Carlo studies showed that the energy resolution of the 
calorimeter would be degraded to 60-70%/\/E'- This large sampling term will 
dominate the energy resolution even at the HERA kinematical end point of 
820 GeV provided the constant term is held to 2% or better. 

2.2 Description of the Calorimeter 

The FNC, built in two parts (Figs. and occupies a very restricted space 
in the HERA tunnel between the electron beam line and the tunnel wall, after 
the 6 mrad upward bend of the protons so that they can circulate above the 
electrons. It lies immediately in front of the beginning of the cold section of 
the proton machine. The available space beneath the proton beam pipe is 
limited in height to 50 cm. The zero degree line enters the calorimeter 16 
cm below the beam pipe. The top part of height 20 cm was constructed to 
cap the bottom part by surrounding the beam pipe. The two parts together 
provide symmetrical up-down coverage, 35 cm in diameter. The total length, 
including readout and scintillation counters, is limited to 3 m. Since space is 
needed for PMT assemblies and scintillation counters both before and after 
the calorimeter, only about 2 m was available for the active calorimeter itself. 

The calorimeter is 134 layers deep with a unit sampling cell as given in 
Table 2. The energy sampling fraction is 3% for minimum ionizing particles 
(mip) and 2% for electrons, as determined using GEANT. For readout, it 
is divided longitudinally into a front section of 95 layers and a rear ‘tail 
catcher’ of 39 layers. This division serves two purposes. By shortening 
the length of the wavelength shifter, it reduces the difficulty of ensuring 
the longitudinal uniformity of response (section |2.3|) . By having a shower 
tail catcher, late starting showers which have increased probability of large 
longitudinal energy leakage can be tagged. Since the rear is not needed for 
position measurements, its transverse segmentation is reduced. 

The calorimeter is read out on both sides by wavelength shifting light 
guides (WLS) and photomultiplier tubes (PMT). The sections of the calorime¬ 
ter, namely, bottom front and rear, and top front and rear, have different 
widths. The maximum width of 74 cm was constrained by the presence of 
electron quadrupoles and radio-frequency cavities. Table ^ gives the nominal 
size of each section of the ENG. 
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Material 

Depth 

(cm) 

Absorption Lengths 
Aa 

Pb 

1.25 

0.073 

scintillator 

0.26 

0.003 

paper 

0.04 

0 

air 

0.1 

0 

total 

1.65 

0.076 

134 layers 

221.1 

10.2 


Table 1; The structure of a unit cell of the FNC. 


Section 

layers 

length 

(cm) 

depth 

Aa 

height X width 
(cm X cm) 

tower size 
(cm) 

PMT channels 

Bottom front 

95 

156.75 

7.2 

50 X 70 

5 

20 

Bottom rear 

39 

64.35 

3.0 

50 X 74 

10 

10 

Top front 

95 

156.75 

7.2 

20 X 60.5 

5 

8 

Top rear 

39 

64.35 

3.0 

20 X 64.5 

10 

4 


Table 2: The size of the FNC sections. All the scintillator strips have a true 
height of 4.94 cm. 
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The calorimeter readout is divided into vertical towers to permit a mea¬ 
surement of the impact position of the scattered neutron on the face of the 
calorimeter. This position is needed for the measurement of pt of the scat¬ 
tered neutron. The vertical position, y, is determined by energy sharing 
between towers; the horizontal position, x, by light (collected charge) divi¬ 
sion between the left and right sides; the HERA proton beam has an average 
angular spread at the interaction point of about 70 prad. This results in 
an intrinsic spread in the neutron impact position of 0.7 cm, independent of 
its Pt- This intrinsic spread sets the limit for position resolution. With the 
assumption that the position resolution scales as 1/V^, and the requirement 
that neutrons with E> 200 GeV have a position resolution better than 0.7 
cm, we demand a resolution of 10 cm/v^- Experience with the uranium 
calorimeter [pTjl in ZEUS indicates that this can be achieved, vertically, with 
tower segmentation of 5 cm. The scintillator, SCSN-38, was cut into strips 
4.94 cm high, and wrapped in white Tyvek paper. All the scintillator strips 
have the same height, including the ones in the rear section where vertical 
segmentation of the readout is 10 cm. Here each wavelength shifting light 
guide reads out two rows of scintillator strips. 

The horizontal position resolution depends on n^e, the number of photo¬ 
electrons collected, and on A, the attenuation length of light in the scintillator 
strips. There is a contribution to the position resolution from shower fluctu¬ 
ations. For electrons, this is much less than the photostatistics contribution, 
while for hadrons we expect the fluctuation component to be larger, but the 
photostatistics component still dominant. The error on the horizontal posi¬ 
tion measurement is then 5x = X/ jUveW^- The ZEUS uranium calorimeter 
collects 220 pe/GeV, and ZEUS testsf^ suggest that A would be about 100 
cm for the 70 cm x 5 cm x 0.26 cm scintillator strips. By assuming that 
the light collecting components of the ENG have similar optical properties to 
those of the ZEUS calorimeter, we estimate Upe for the ENG from the ZEUS 
calorimeter number. Since the ENG has a smaller sampling fraction (3:9), 
longer tower widths (70:20), and transmission rather than backing masks 
for the WLS, we estimated the ENG would collect only 45 pe/GeV, for a 
resolution in x of about 15 cm/\/E. 
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2.3 Wavelength Shifter and Photomultipliers 

An excellent wavelength shifting light guide material to match the SCSN-38 
scintillator is Y7 in PMMA base (30 ppm doping); however, the manufac¬ 
turer, KurarayQ had only enough stock on hand for the rear section of the 
calorimeter, so the front section light guides were produced from similar ma¬ 
terial (BC-482A with 20 ppm doping) manufactured by Bicronj^. 

The WLS in FNC were made 0.4 cm thick to absorb as much as possible 
the light leaving the scintillator strips while still having the light guides match 
the size of the PMTs. So that longitudinal shower fluctuations do not degrade 
the energy response of the calorimeter, transmission masks were positioned 
between the WLS and the scintillator to correct for longitudinal attenuation 
of light along the WLS. The masks were a pattern of black lines printed on 
3M overhead projector transparencies using a 300 dpi laser printer. 

The attenuation lengths of the light guides were measured using equip¬ 
ment and procedures developed for the ZEUS calorimeter [^. Figure ^ shows 
a comparison of the curves, before and after masking, for a sample light 
guide. For the front WLS the masking reduces the peak light yield at the 
PMT by 2.2; for the rear, the reduction is by 1.4. The transparency material 
itself accounts for some of this reduction since, as measured by a spectropho¬ 
tometer, it transmits only 79% of incident light at a wavelength of 420 nm 
where the emission spectrum of SCSN-38 peaks. 

After masking the response should be flat. A measure of the residual 
non-uniformity is the root mean square (rms) deviation from a constant of 
the corrected light curves. This averaged 2% of the mean. 

Two types of photomultiplier tube were used: Phillips XP 2008/UB and 
Valvo XP 2011/03. Fifty pairs of tubes and resistive bases, which were spare 
from the ZEUS calorimeter tests at CERN, were obtained and tested in the 
ZEUS photomultiplier tube test facility to determine their characteristic gain 
and linearity. 


2.4 Source Scans 


We studied the FNC by scanning it with a 1 mCi ®°Co source, a procedure 
developed by the ZEUS calorimeter group [T^. This allowed us to check the 


'^Kurary Co., Tokyo, Japan. 

^Bicron Corporation, Newbury Ohio, USA. 
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quality of the transmission masks, to look for construction problems, and to 
measure the relative gain of each channel. 

A mechanical assembly equipped with a stepping motor moved the source 
through brass guide tubes hxed outside the calorimeter along the WLS. An 
IBM 486 PC containing an integrating ADC card controlled the stepping 
motors and measured current from the two PMTs viewing the tower along 
which the source was moving. Readings were taken about 0.15 cm apart. 

Figure ^ shows a typical scan of a front tower, displaying the ADC out¬ 
put as a function of longitudinal position for the PMT on the side of the 
calorimeter opposite the source. Each scintillator layer is visible as a peak. 
As the source is withdrawn from the calorimeter, the response falls and flat¬ 
tens to become the pedestal which must be subtracted to obtain the gain. 
The steepness of the falloff at the edge of the calorimeter and the depth of the 
valley between peaks indicate that the general response is a superposition of 
several scintillator layers, as expected from the low energy (1.2 MeV) gamma 
rays emitted by ®°Co. 

To measure the pedestal and scintillation signals, a scan is projected 
onto the ordinate of the plot. The relative gain of a PMT, compared to all 
others, is obtained by comparing the htted average peak position for its scan 
after subtraction of the pedestal. After correction for pedestal and relative 
gain, the projection of the peaks for all the scintillator layers in all towers is 
distributed approximately as a Gaussian, whose width to mean ratio (4.7%) 
is a measure of the uniformity of transmission of the optical system, including 
scintillator and masked WLS. 


3 Test of the Calorimeter at CERN 

Only the bottom part of the calorimeter was ready for testing at CERN in 
the summer of 1994. The calorimeter was installed in the H6 beam in the 
North Area of the SPS. This beam provided electrons, muons, pions and 
protons up to an energy of 120 GeV. The positively charged particle beam 
is mixed pions and protons with no particle identihcation[|T^. Most of the 
data was taken at 120 GeV. 

The calorimeter was placed on a table whose horizontal and vertical po¬ 
sition were adjustable and measurable to ±0.1 cm. A pair of crossed huger 
counters in front of the calorimeter was used for triggering and to dehne the 


location of the beam (Table Three large counters were also positioned in 
front of the calorimeter. The signal from the last dynode of each of the eight 
PMTs of the four central calorimeter towers was also available for triggering 
on calorimeter energy deposits. 


Counter 

Description 

Coverage 

T1 

vertical huger counter 

1 cm X 0.5 cm 

T2 

horizontal huger counter 

1 cm X 0.5 cm 

T3 


20 cm X 20 cm 

T4 


15 cm X 15 cm x 5 cm 

T5 

muon counter 

20 cm X 20 cm 

T6 

veto 

3 cm diameter hole 


Table 3: Beam counters. 

The photomultipliers were read out with LeCroy 4300 charge integrating 
ADCs. The ADCs were gated so that 160 ns of the pulse from the calorime¬ 
ter was integrated. During data taking the ADC pedestals were regularly 
measured, and the gains were monitored with an internal test function in the 
units. The gains were also checked daily using an external charge injector. 
Over the course of the beam test, which lasted one week, individual ADC 
gains varied less than 1 part per mil, and all gains were uniform to within 5 
parts per mil. 

3.1 Coordinate System 

In the following we use a right-handed coordinate system with y vertically 
up and 2 ; into the face of the calorimeter. The origin (0,0) is taken to be 
between the two central towers halfway between the sides (the center of the 
face of the calorimeter). X and Y are measurements of x and y based on the 
table position; Xcal and Tcal, measurements based on energy deposits in 
the FNC. 
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3.2 Measurement of Shower Position and Size 


The vertical position of a shower is determined using the centroid method 
by taking an energy weighted average of the the tower positions: 


WiVi 


dcAL = 




( 1 ) 


where the sum is over towers in the front part of the FNC, and the weights Wi 
are functions of the energy deposits, Ei, in each tower. The weights are chosen 
to minimize bias and to maximize resolution in the position measurement. 
The simplest method of estimating y corresponds to choosing weights which 
are linear in energy, Wi = Ei. Awes et al.|^ have suggested that the weights 
be chosen to have instead a logarithmic dependence on the energy deposits: 


Wi = max 



( 2 ) 


where Eb = Ei, and 0 < / < 1 is a fractional cutoff. When the 
energy deposited in a tower satishes E^ < fE-Q, that tower is not used in the 
determination of the vertical position. The logarithmic dependence of the 
Wi more closely matches the true transverse prohle of energy deposits in a 
hadronic shower than the linear dependence. This estimator is simple, but 
it has an undetermined free parameter / which may be a function of the 
deposited energy. In addition the optimal value of the cutoff parameter may 
depend on particle type, electron or hadron. 

The vertical segmentation also allows us to estimate the size of the shower 
using the RMS vertical shower width: 


W 


E,‘£l 


\ EZ, 


Wi 


^CAL- 


(3) 


Since the width of electromagnetic showers is much less than the width of 
hadronic showers, this estimator enables us to distinguish incident electrons 
from incident hadrons. The optimal weights for the width determination are 
not necessarily the same as those for the position measurement. Since we 
use the width W to separate electromagnetic from hadronic showers, it is 
most convenient to have an estimator that does not depend on knowledge of 
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the nature of the incident particle. Otherwise an iterative procedure must 
be employed. As we will discuss in more detail below, linear weights for the 
shower size measurement give a hadron-electron separation which is close 
to that obtained with logarithmic weights whose cutoff parameter has been 
optimized. 

The horizontal position of a shower can be estimated from the division of 
scintillation light between the left and right sides of the calorimeter: 

= (4) 

where A is the effective attenuation length of the scintillator strips, and L{R) 
is the charge from the left (right) PMT. If the light attenuation behaves as 

L,R (X exp{±ax + bx‘^), (5) 

where a, b are constants, and x is the displacement of the incident particle 
from the center of the tower, then Xcal is independent of b and A = 1/a. 

3.3 Electrons 

The calorimeter was calibrated with 120 GeV electrons incident on the center 
of each tower. The high voltage for the tubes in the front was adjusted so 
that their mean response to 120 GeV electrons as measured in the Lecroy 
ADGs was 100 pG (400 counts). After lifting and rotating the calorimeter by 
180 degrees to make electrons incident on the rear towers, the high voltage 
of each rear tube was set so that the response of the rear section to a muon 
passing through the calorimeter would approximately equal the response of 
the front section. Since the ratio of the number of front layers to the number 
of rear layers is 95/39, the gain is considerably higher in the rear; however, 
the energy deposits in the rear for electrons and hadrons are small, so there is 
no saturation problem for the very high energy particles at HERA. Figure 
shows the distribution of deposited energy, after calibration, for 120 GeV 
electrons centered on towers 3 through 7, inclusive. Pion contamination has 
been removed be requiring that the vertical shower width be less than 3 cm 
(section |3]^). The energy resolution for electrons centered on a tower aver¬ 
aged 33%/independent of tower. Electrons centered in the edge towers 
had an energy response reduced by 3% due to transverse energy leakage. 
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The response of the calorimeter to electrons is dependent on the posi¬ 
tion of the electron with respect to the gap between the scintillator strips. 
Fignre ^a shows the energy deposited in tower 6 when the electron beam 
is scanned from tower 6 to tower 7. The total signal for electrons incident 
between two towers is 90% of the signal of an electron incident on the center 
of a tower. Moreover, the energy resolntion is degraded to 60%/a/E'. These 
effects are dne to the narrowness of electromagnetic showers compared to the 
width of the region between strips. 

The small width of the shower compared with the strip width, combined 
with the rednction in response at the tower bonndaries, is responsible for 
the large bias towards the tower center in the measnred electron vertical 
position shown in Fig. The bias grows with increasing cntoff parameter 
since for snfficiently large / only the strnck tower contribntes to the snms 
in eqnations and |^. The RMS deviation of the distribntion also tends 
to increase. The variation with vertical position of the shower width (linear 
weights) of electrons is shown in Fig. and is small compared to the average 
width of hadronic showers (see Section |3.5| ). 

Electrons were also used to determine the effective attenuation length for 
light in the scintillator strips. One half of each tower was scanned horizon¬ 
tally. For each position the response of the right and the left photomultiplier 
tubes was plotted as a function of the offset x of the beam from the center of 
the calorimeter, positive x representing an offset away from the PMT. Fig¬ 
ure H shows the resulting attenuation curve for a typical tower. The response 
is well represented by equation |^. The tower-to-tower average value of a was 


found to be 1.08T0“^ cm“^ with a 3% RMS spread, corresponding to an av¬ 
erage effective attenuation length of 92 cm; b averaged 9.60T0“^ cm“^ with 
a 20% RMS spread. 

To compare the efficiencies of the wavelength shifting material in the front 
and rear section we measured the number of photoelectrons collected by the 
photomultiplier tubes, ripe, which is related to the width, a{L — R), of the 
difference distribution by 


a{L -R) _ 1 

< L + R> VUpe ’ 


( 6 ) 


where < > is the mean of the sum. By rotating the calorimeter 180 de¬ 

grees to make electrons incident both front and rear, we found that although 
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in the rear we averaged 92 pe/GeV, in the front we obtained only 26 pe/GeV. 
This difference is due to different WLS lengths, WLS material, and optical 
geometry between front and rear, and it worsens the x position resolution 
by a factor of a/2 compared to the estimate in section p.2| . Photostatistics 
makes only a small contribution to the energy resolution since, for 120 GeV 
electrons and 26 pe/GeV, the fluctuations are only about 2%. 

3.4 Hadrons 

Electron showers penetrate only about an interaction length into the calorime¬ 
ter. To rehne the calibration constants determined by electrons, 120 GeV 
pions were made incident on the center of each tower. Using these data, cor¬ 
rection factors were determined which equalized the response of each PMT 
to hadrons. The correction factors for hadrons were forced to average 1.0, 
and had an RMS spread of 5%. The distribution of deposited energy for 120 
GeV hadrons (pions and protons) incident on the center of the calorimeter, 
between towers 5 and 6, is shown in Fig. For pions incident in the center 
of each tower, the average response of the calorimeter is shown in Fig. ^ The 
curves are quadratic in y^. Leakage loss is important whenever the beam is 
outside the four central towers. When the beam is incident in the center of 
the calorimeter, between towers 5 and 6, about 2.5% of the energy is lost 
compared to the case when the beam is in the center of tower 5 or 6 because 
of the gap between the scintillator strips (see inset of Fig. |^. Parameterized 
as ct/a/E, the energy resolution for hadrons incident at the center of towers 
5 and 6 is 0.70/V^ (Fig. This improves to 0.62/a/E for hadrons incident 
at the center of the calorimeter, between towers 5 and 6. 

These results can be compared to the results for electrons discussed above. 
The electron to hadron ratio e/h is 0.96 when the beams are incident on the 
center of a tower, but it falls to about 0.9 when both are incident on the center 
of the calorimeter in the gap between towers 5 and 6. For electrons centered 
on a tower, but hadrons incident on the center of the calorimeter e/h = 0.99 
(see Fig. |). If the electron calibration constants are used instead of the 
hadron calibration constants, e/h changes by 0.01. As expected from the 
ratio of lead to scintillator thickness, the calorimeter is over compensating. 
This effect is magnihed for electrons by the gap between the scintillator strips; 
however, since compared to the electrons the electromagnetic component of 
a hadronic shower is more spread out transversely, we expect the ‘effective’ 
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e/hiohe higher. 

In order to study energy leakage and position resolution for hadrons, 
the face of the calorimeter was scanned with 120 GeV pion over an x, y 
grid. Using this grid the effective attenuation length of the scintillator was 
determined with hadronic showers. The resulting average attenuation length, 
91.9 cm, was in good agreement with the electron scan results. 

Figure shows the true vertical position of the calorimeter as a function 
of its position determined using the centroid algorithm with linear weights 
{wi = Ei). The effect of the hnite strip size is clearly evident in the os¬ 
cillation superimposed on an average behavior. In addition, towards the 
edge of the calorimeter the measured value falls below the true value be¬ 
cause of leakage. The average behavior can be parameterized conveniently 
as |U| = I.OSOITcalI + 0.002y3AL- The residuals, shown in the upper inset, 
follow an approximate sine curve with period 5 cm, the tower height. 

The overall vertical resolution is determined by plotting the y residuals 
for the centroid algorithm with the polynomial correction for all points of the 
grid. The residual distribution is well £t by a sum of two Gaussian (Fig. 0a). 
The largest component has a weight of 0.91 and an RMS width of 0.8 cm for 
120 GeV pions. Because this width is larger than the size of the oscillation, 
the y position cannot be corrected for the oscillation on an event by event 
basis. 

In contrast to linear weights, logarithmic weights depend on an arbitrary 
cutoff parameter / which may be chosen to optimize the position resolution, 
minimizing bias. Residuals for logarithmic weights with / = 10% are shown 
in the lower inset in Fig. 0 and in Fig. 0b. The bias due to leakage has 
been considerably reduced (|V| = O.OOlycALl + O.OOlVcAL), and the resolution 
improved, compared with the use of linear weights. 

A measure of the effective overall resolution is (Teff 

‘^EFF ~ + f2(^2 ('^) 

where af is the variance of Gaussian z, and fi is its fractional contribution to 
the overall distribution. The resulting cteff is plotted in Fig. 0 as a function 
of /. The best vertical position resolution for 120 GeV pions is 8 cm/ a/U for 
/ = 9.4%, where the curve has a minimum. 

Like the y residuals, the distribution of x residuals is also well represented 
by a sum of two Gaussians: 97% with cr = 20 cm/ a/E, and 3% with cr = 60 
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cm/V^, giving cteff = 22.3 chi/a/E. The resolution in x is much poorer 
than in y, nevertheless, it is close to the value expected since for 1 GeV 
incident energy 


6x = \j^ 92 cm./\f{2Q) = 18 cm. (8) 

The energy response varies over the grid as a result of two competing 
effects: a) energy leakage, and b) light enhancement. As the beam is moved 
towards the edge of the calorimeter, energy is lost due to leakage, and the 
measured signal tends to decrease. On the other hand, as the beam moves 
horizontally away from the center the light collected by the PMTs increases 
because of the characteristic light attenuation of the scintillator strips: 

L + i? oc exp(6a;^) (exp(—ax) + exp(ax)). (9) 

The net result is that the signal increases for horizontal displacements at 
hxed y. The increase is approximately linear in x with 

L + R(x[l + (0.002 + 0.00006||/|)|x|]. (10) 

Figure |T^ shows the energy response as a function of y for all grid values 
with the X dependence removed. Energy loss due to the gap between the 
scintillator strips is clearly visible as is also loss due to leakage. The leakage 
component has a quadratic dependence in y'^; the loss due to the gaps between 
strips behaves as a superposed oscillation. 

3.5 Electron-Hadron Separation 

A calorimetric measurement with the capability of distinguishing between 
electromagnetic and hadronic showers is important for physics studies. Be¬ 
cause electromagnetic showers are much narrower than hadronic showers, the 
measurement of the shower width allows a separation of particles by type. 
The use of logarithmic weights for the determination of the vertical shower 
position reduces the bias due to the hnite strip size and lateral leakage, and 
improves the resolution; however, the optimal choice of the cutoff parameter 
depends on incident particle type. At 120 GeV, / should be chosen to be 
about 10% for hadrons, but should be close to 0 (no cutoff) for electrons. 
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The choice of cutoff for the width measurement can be chosen to optimize 
the separation of electromagnetic and hadronic showers. To estimate the 
optimal value of / we define 


Sep = 


< Wh > - < We > 


( 11 ) 


as a measure of the electron-hadron separation, where < Wh,e > and Sh,e are 
the means and standard deviations, respectively, of the width distributions 
for hadrons {h) and electrons (e). All are implicitly functions of the cutoff 
/, so the value of / can be chosen to maximize Sep. 

The dependence of Sep on / for 120 GeV electrons and hadrons is shown 


in Fig. |TJ. Sep is insensitive to / at low values of /, but falls approximately 
linearly for />5%. Linear weights give the result shown by the dashed line, 
Sep = 3.4, close to best values obtained with logarithmic weights, and have 
the advantage of being independent of /, which may depend on the incident 
energy. As a result, we use linear weights for measurement of shower widths 
instead of the logarithmic weights chosen for the shower centroid measure¬ 
ment. 

The vertical width distribution is plotted in Fig. 0b and c for 70 and 
120 GeV electrons and pions incident on towers 2 through 9. Although the 
distributions do not change much with energy, there is a small reduction 
in spread and decrease in mean. Figures igc and p!5| a show that mean and 
standard deviation of the shower width distribution for hadrons does not 
change significantly with the incident vertical position. A W cut chosen in 
the range 2-4 cm will clearly separate electrons from hadrons. Only 2.3% of 
hadrons incident at the center of the calorimeter will be lost with a W cut 
of 3.5 cm; the loss rises to 4.0% for hadrons incident at the center of a strip 

(Fig. 0). 

The same analysis for electrons is complicated by pion contamination in 
the electron sample, which is shown as a broken line in Fig. 0 To subtract 
the pion contamination, we normalize the width distribution for pions to the 
raw (unsubtracted) distribution for electrons in the region W > 6 cm. The 
result is shown in Fig. p!5[ b for electrons incident in the center of a tower, and 
electrons scanned vertically from tower 6 to 7. The inset shows the tail of 
the distribution for centered electrons. About 1.5% of the electron showers 
have a width greater than 3.5. 
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3.6 Muons 


Although the signal from penetrating muons is low because of their small 
energy deposits, they can be useful for calibrating and monitoring the FNC 
in HERA using either cosmic rays or beam halo muons. To study the response 
of the calorimeter to high energy muons, a special trigger was implemented for 
beam particles which penetrated the calorimeter. Figure IB shows a typical 
energy distribution for a muon in the bottom front tower. The peak position 
Xo is determined by fitting the distribution to a convenient functional form: 


m{x) = Aexp (a (?/— exp(—?/))), 

X — Xq 


( 12 ) 


V = 


a 


and A, a, xq, and a are parameters. When a = 1/2, m{x) is the Moyal 


function|T^, an analytic approximation for the Landau distribution. The 


tower-to-tower spread of the peak position is 4%. 


3.7 Neutrons 

The absolute energy scale is best determined by making a beam of known 
energy incident on the calorimeter. If deuterons were to be accelerated at 
HERA, nuclear stripping reactions would make neutrons of known energy 
incident on the calorimeter. At present only protons are accelerated, and the 
absolute energy scale must be determined using beam gas interactions which 
produce leading neutrons by charge exchange. If deuterons were accelerated 
at HERA in order to test the feasibility of using charge exchange reactions 
to determine the absolute energy scale, a short experiment was performed. 

Figure ^ shows the experimental arrangement. A Incite target was placed 
in the beam 20 m in front of the calorimeter. Immediately after the target a 
sweeping magnet was tuned to deflect beam particles downward. With the 
magnet tuned to full field beam particles ( protons and pions ) were deflected 
downwards by 11 cm at the calorimeter. It was necessary to devise a trigger 
that ensured that only events arising from the impact of a neutral particle 
in the central part of the calorimeter were accepted. 

Crossed finger counters in front of the target defined the beam position 
(see Table and Fig. 0)- The trigger required a coincidence of these two 
counters and energy deposited in the center of the calorimeter For the trig¬ 
ger, the deposited energy was measured by summing and discriminating the 
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signals from the last dynode of the PMTs of the four central towers. In order 
to veto charged particles, three scintillation counters were placed directly in 
front of the calorimeter. In addition, only showers near the center of the 
calorimeter were considered in order to reduce energy leakage. In summary 
the data were required to satisfy: 


(T1 > mip) ■ (T2 > mip) (13) 

(T3 < mip) ■ (T4 < mip) ■ (T6 < mip) (14) 

Si^dynodes > threshold (15) 

R < 5 cm, (16) 

where 

R = Aoal + YSal- (17) 


The shower position Xcal was determined by light division; Ycal, by the 
centroid method with logarithmic weights. 

The observed energy distribution is shown in Fig. |T^. In the same hgure, 
the hatched histogram shows the response of the calorimeter to 120 GeV 
hadrons. The mean value of the latter distribution determines the energy 
scale for the hgure. 

The one pion exchange (OPE) model of the reaction p ^ n predicts that 
the cross section, integrated over scattering angles, is given by 


da 

dE 




(18) 


where /^r/p is the hux of virtual pions associated with the incoming proton; 
cjTr is the total pion-target cross section for a virtual pion beam energy of 
AE = Efj — E. Any dependence of on the virtuality of the exchanged pion 
has been ignored in equation so a.,, can be estimated using a parameteri¬ 
zation {AE > 4 GeV), or an interpolation table {AE < 4 GeV)|^. Several 
forms have been postulated for the hux factor /^/p in theoretical and phe¬ 
nomenological studies of OPE. They fall into two classes: Regge[p!^, 
and light cone . 


Accordingly, to compare data and theory, we extended the IIERWIG| 
Monte Garlo program to generate events according to OPE folded with the 
measured experimental resolutions. The predicted energy distribution, for 
the light cone hux factor]^, adjusted for resolution and acceptance is shown 
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as the open histogram in Fig. |^. The OPE prediction is normalized to the 
data by the number of events. As can be seen, agreement is good for energies 
larger than 100 GeV. 

The energy scale was varied and then compared with the OPE prediction 
for E > 100 GeV by computing the of fh® difference distribution. The 
variation of with the scale factor is shown in the inset to Fig. |^. The 
minimum occurs for a scale factor of 0.985. If the Regge |T^, ^ form for the 
flux factor is used instead, the same procedure yields a minimum at 1.003, in 
complete agreement with the incident proton beam; in this case, however, the 
minimum is worse (17 rather than 10). We conclude that using charge 
exchange neutrons the absolute energy scale can be determined to better 
than 1.5% with the chief source of error being the theoretical uncertainty on 
the form of the virtual pion flux factor. 


4 Monte Carlo Studies of the FNC 

Since the spring of 1995, the FNG has been operating in the ZEUS experiment 
at HERA where it is used in the study of leading neutrons produced at 
small angles with energies E> 100 GeV. Because the calorimeter has been 
calibrated and tested only in beams of energy up to 120 GeV, and because 
the top was not present for the beam tests, we must rely on a Monte Garlo 
simulation to predict the response of the calorimeter to high energy particles. 
We have modeled the FNG using the GEANT 3.13|Q program, upon which 
the simulation of the full ZEUS detector is based. In this section we present 
some results from the simulation which can be compared to our test beam 
data. 

For 120 GeV electrons and pions incident on the center of towers 5 and 6 
the GEANT simulation predicts an electron to hadron response ratio of 0.98, 
in agreement with the measured value of 0.96. The simulated response to 
pions incident on the center of each tower is shown in Fig. |T^. The energy loss 
due to leakage when the beam is incident near the edge of the calorimeter is 
also in agreement with the data as is the degradation of the energy resolution 
(compare with the data shown in Fig. P). 

The Monte Garlo gives an energy resolution due to shower fluctuations 
alone, that is, neglecting fluctuations due to photostatistics, of 0.66/-\/U for 
hadrons incident at the center of towers 5 and 6. We have also used the 
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Monte Carlo to predict, at higher energies, the overall energy response of the 
calorimeter and its expected energy resolution due to shower fluctuations. 
The Monte Carlo predicts that for neutrons centered on the calorimeter with¬ 
out the top, the energy response is linear up 800 GeV. Fitting the energy 
resolution as a function of incident energy, we find 




0.54 

Ve 


© 0.03 


(19) 


where E is in GeV. If fluctuations due to photostatistics are added, the 
sampling constant 0.54 will increase to 0.58. The top section is present at 
HERA, but the neutrons are predominantly incident on towers 7 and 8. 

Just as for data, 120 GeV pions were studied with a grid over the face 
of the calorimeter. Fig. ^ shows the predicted dependence of Tcal on V, 
the vertical impact position, for three values of horizontal impact position, 
X = 0, 10.1, and 20.1 cm. Ycal is calculated with logarithmic weights and 
a cutoff parameter of f=10%. The simulation shows behavior and biases 
similar to the data (Fig. 0 lower inset); in particular, the value of Tcal 
is biased towards the nearest tower center. The Y residual distribution is 
Gaussian with mean 0 and width 7.3 cm./\fE (compare to the data width 
of 8.0 cm/\fE). The X residual distribution is Gaussian with mean 0 and 
width 10.3 cm/\fE. This width is due to transverse shower fluctuations. If 
photostatistics fluctuations are also included, the width constant increases 
to 20.7 cm. This should be compared with the data value of 22.3 cm. 

The simulated energy response for pions as a function of position, shown 
in Fig. is in agreement with the data shown in Fig. Fig shows 
simulations of the vertical shower widths for 70 and 120 GeV electrons and 
pions. The Monte Garlo results are narrower than the data shown in Fig. |T^ 


5 Summary 

We have designed and constructed a lead scintillator sandwich calorimeter for 
the ZEUS experiment at HERA. The calorimeter is divided into 5 cm vertical 
towers read out on two sides with wavelength shifting light guides coupled 
to photomultiplier tubes. The calorimeter was tested in the H6 beam line 
at GERN with 120 GeV electrons, muons, pions and protons. Electrons can 
be cleanly separated from hadrons using the energy weighted vertical width 
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of a shower. At 120 GeV the calorimeter is slightly over compensating, with 
an electron to hadron response ratio of 0.96 and an energy resolntion 6% 
at 120 GeV. The horizontal position resolntion, measnred by charge sharing 
between the two sides, is 20 cm/v^, assnming a 1/Ve dependence; the 
vertical position resolntion, measnred by energy sharing between the towers, 
is 10 cm/y/E. By nsing energetic nentrons prodnced by proton interactions 
in a Incite target, the overall energy scale was determined to 1.5%. 

Since the spring of 1995, the calorimeter has operated snccessfnlly in the 
HERA tunnel 105.6 m downstream of the ZEUS detector on the zero degree 
line. 
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Figure 1: An isometric view of the FNC showing the top and bottom front 
and rear sections with their photomultiplier tube assemblies. At HERA 
the neutron line of flight, from the bottom left to the top right, enters the 
calorimeter 16 cm below the top of the bottom section. The top section was 
not present for the CERN tests. 
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Figure 2: A schematic view of the front of the FNC. For the CERN test, the 
front towers in the bottom part were labeled from 1 to 10, starting at the 
bottom. All dimensions are in mm. 


25 











Co60 Scan of the FNC 



Figure 4: A far side ®°Co source scan of tower 3 of the front section of the 
FNC. The peaks of the 95 scintillator layers are clearly visible. The inset 
shows an enlargement of a short section. 
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Events Events 




Figure 5: The energy distribution in ADC counts for (a) 120 GeV electrons 
centered on towers 3 to 7, and (b) 120 GeV pions centered on the calorimeter, 
between tower 5 and 6. Gaussians are fitted to the distributions. The pion 
contamination in the electron sample has been removed by requiring W < 3 
cm (see text). 
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Figure 6: (a) The average energy, in ADC counts, deposited in in the central 
towers, 6 and 7, as a function of Y position for a vertical scan from tower 6 
to tower 7. The point at Y = 2.5 cm, the center of tower 6, is hxed at 800 by 
the calibration with electrons. The gap between towers is near Y = 5 cm. (b) 
The average width measured for 120 GeV electrons in a vertical scan across 
the gap between the scintillator strips of tower 6 and 7. The data points are 
means; the error bars, the corresponding RMS deviation of the distributions. 
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120 GeV Electrons 
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Figure 7: The average measured Tcal as a function of Y for 120 GeV elec¬ 
trons in a vertical scan across the gap between scintillator strips. Shown 
is VcAL as measured using (a) linear weights and using logarithmic weights 
with a cutoff of (b) 10% and (b) 1%. The data points are means; the error 
bars the corresponding RMS deviation of the distributions. 
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120 GeV Pions 



Figure 10: y as a function of Fcal, computed with linear weights, for 120 
GeV pions incident on the face of the calorimeter over a grid in x and y (see 
Fig. |T^ for the grid and dehnition of symbols). The upper inset shows the 
residuals from a quadratic fit. The lower inset shows the residuals when Fcal 
is computed with logarithmic weights (f=10%). The dashed line is V = Fcal- 
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120 GeV Pions 




(b) Y-Yc^ (cm) 


Figure 11: The distribution of residuals of the measured y position for 120 
GeV pions incident on face of the calorimeter a grid in x and y (see Fig. for 
(a) linear weights, and (b) logarithmic weights (f=10%). The distributions 
are fit to the sum of two Gaussians. 
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Figure 12: The variation of the effective vertical position resolution, deter¬ 
mined with logarithmic weights, as a function of the cutoff fraction /. The 
minimum of the curve, where the resolution is best, occurs at / = 9.4%. 
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Electron-Hadron Separation 
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(c) Width (cm) 

Figure 14: (a) The electron-hadron separation variable, Sep, as a function 
of the logarithmic weight cutoff parameter /. The dashed line shows the 
value of Sep for linear weights. Below (a), the distribution of shower widths, 
computed with linear weights, for (b) 70 GeV and (c) 120 GeV electrons 
and pions. The electrons are incident on the center of towers. The pions are 
incident at the center of a tower (full) and at the center of the calorimeter 
(dotted). The pion contamination in the electron sample is visible. 
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Figure 15: (a) The variation of the mean and RMS of the width distribution 
of 120 GeV pions with vertical and horizontal position in the calorimeter. 
The data points are the mean of the distribution; the ‘error’ bars the RMS. 
(b) The distribution of shower widths, computed with linear weights and 
after subtraction of the pion contamination, for 120 GeV electrons centered 
on a strip (full), and a scan of electrons starting from 0.9 cm above the center 
of tower 6 to tower 7 (dashed, see Figs. ^ and |^). The inset shows the tail of 
the distribution for centered electrons. 
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Figure 17: A schematic of the FNC test setup at CERN. 




120 GeV p^n 



Figure 18: The neutron energy spectrum (solid points) compared to the 
spectrum predicted by one pion exchange (open histogram). The hatched 
histogram shows the measured energy spectrum for 120 GeV hadrons (see 
Fig. |b). The inset shows the comparison of the measured neutron spec¬ 
trum and the predicted neutron spectrum as a function of the neutron to 
beam hadron energy scale factor (sf). 
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Figure 19: Monte Carlo simulations of the response of the calorimeter to 120 
GeV pions incident on the center of each tower. The energy resolution, due 
to shower fluctuations, that is, without photostatistics, is shown in the lower 
graph. The corresponding test beam data are shown in Fig. p. 




0.667/VE 


I 


I 


♦ * ♦ ♦ r, 




I I 1 1 I I 1 I I 1 1 1 1 I I 1 X I 1 1 I 


42 





120 GeV Pions (MC) 



Figure 20: Monte Carlo simulation of Fcal, determined by logarithmic 
weights with f=10%, as a function of Y for 120 GeV pions incident on the 
calorimeter over a grid in x and y. See Fig. 0 for the corresponding plot 
made with test beam data. 
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E (Counts) 


120 GeV Pions (MC) 



Figure 21: Monte Carlo simulation of the measured energy as a function of 
Y for 120 GeV pions incident on the face of the calorimeter. See Fig. for 
the corresponding plot made with test beam data. 
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Figure 22: The histograms are Monte Carlo simulations of the width distri¬ 
butions, with linear weights, for (a) 70 and (b) 120 GeV electrons and pions 
incident on tower 6. See Fig. for the corresponding plot made with test 
beam data. 
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